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^ ■ We compare updated predictions for the elastic scattering of supersymmetric neutralino 

p • dark matter with the improved experimental upper limit recently published by CDMS II. We 
take into account the possibility that the 7r-nucleon S term may be somewhat larger than 
! was previously considered plausible, as may be supported by the masses of exotic baryons 
>• ■ reported recently. We also incorporate the new central value of rrit, which affects indirectly 
^ ■ constraints on the supersymmetric parameter space, for example via calculations of the relic 
I density. Even if a large value of S is assumed, the CDMS II data currently exclude only 
small parts of the parameter space in the constrained MSSM (CMSSM) with universal soft 
supersymmetry-breaking Higgs, squark and slepton masses. None of the previously-proposed 
CMSSM benchmark scenarios is excluded for any value of S, and the CDMS II data do not 
impinge on the domains of the CMSSM parameter space favoured at the 90% confidence 
level in a recent likelihood analysis. However, some models with non-universal Higgs, squark 
and slepton masses and neutralino masses ^ 700 GeV are excluded by the CDMS II data. 
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1 Introduction 



The lightest supersymmetric particle (LSP) is stable in models in which R parity is conserved, 
in which case it is a suitable candidate for the cold dark matter required by astrophysical 
and cosmological observations [1]. One of the generic possibilities is that the LSP is the 
lightest ncutralino x, in which case the detection of dark matter appears feasible. The direct 
detection of supersymmetric dark matter via scattering on nuclei in deep-underground, low- 
background experiments has been discussed many times [2] - [9]. 

There are, however, three reasons why a re-evaluation of the prospects for such ex- 
periments is now timely. The first is the motivation provided by the upper limit on the 
dark-matter scattering cross section provided by the CDMS II experiment [10], which is 
substantially more stringent than previous experiments [11]. The CDMS II result appears, 
in particular, to conflict with the dark-matter scattering interpretation of the results of the 
previous DAMA experiment [12]. A second reason is evolution in Standard Model inputs 
into the calculation of the scattering matrix elements. Recent particle-physics experiments 
tend to favour a value of the pion-nucleon sigma term S that is somewhat higher than ear- 
lier experiments, favouring a larger theoretical estimate for the spin-independent part of the 
dark-matter scattering cross section [13]. Interestingly, a larger value of S is also favoured 
independently by hints from the spectroscopy of pentaquark baryons, if they exist [14]. We 
also include the effect of the new preferred value of rrit [15] on the supersymmetric parameter 
space and on relic-density calculations. Finally, there has been some progress recently in un- 
derstanding which parts of parameter space are favoured in certain versions of the minimal 
supersymmetric extension of the Standard Model (MSSM). In particular, if the input soft 
supersymmetry-breaking parameters are constrained to be universal (CMSSM), the data 
on ruw, sin^ 9w and — 2 all favour independently a relatively low mass for the lightest 
neutralino [16], favouring in turn a relatively large dark-matter scattering cross section. 

It is the purpose of this paper, in light of these developments, to rc-cvaluatc the prospects 
for discovering dark-matter scattering in forthcoming experiments. We include in our anal- 
ysis not only models in which neutralinos are the dominant source of cold dark matter, but 
also those in which neutralinos provide only some fraction < 1. In the latter case, we 
assume that neutralinos constitute the same fraction < 1 of the galactic halo. For com- 
parison with experiments searching for dark-matter scattering, which usually assume that 
all the halos are composed of neutralinos, we rescale the effective scattering cross section by 
the same factor < 1. 

We find that, even with the larger value of E, only very small parts of the CMSSM param- 



1 



eter space are excluded by the current CDMS II result. Specifically, none of the benchmark 
scenarios proposed recently [17] is excluded, and neither is any of the 90% confidence-level 
region favoured in a recent likehhood analysis of the CMSSM [16]. On the other hand, if one 
relaxes universality for the squark slepton and Higgs masses, so as to consider the most gen- 
eral low-energy effective supersymmetric theory (LEEST), some models with ^ 700 GeV 
are excluded for large E. We reach a similar conclusion even if the squark and slepton masses 
are assumed to be equal, and we allow only non-universal Higgs masses (NUHM). Indeed, 
as we discuss, the dominant mechanism leading to a large cross section is the reduction in 
the magnitude of the Higgs superpotential mixing parameter /i and the pseudoscalar Higgs 
mass niA allowed by the relaxed electroweak vacuum conditions in the NUHM. 

2 Spin-Independent x-Nucleon Scattering Matrix Ele- 
ments 

2.1 Model-Dependent Supersymmetric Operator Coefficients 

We assume that the neutralino LSP x is the lightest eigenstate of the mixed Bino B, 
Wino W and Higgsino H12 system, whose mass matrix N is diagonalized by a matrix Z: 
diag{m^^^„^4) = Z*NZ~^, with 

X = Z^iB + Z^2W + Z^sHi + Z^^H2. (1) 

We neglect the possibility of CP violation, and assume universality at the supersymmetric 
GUT scale for the C/(l) and SU{2) gaugino masses: Mi^2 = ?7ii/2, so that Mi = | tan^ 9wM2 
at the electroweak scale. 

The following low-energy effective four-fermion Lagrangian describes spin-independent 
elastic x-nucleon scattering: 

^ = asiXXQiQi, (2) 

which is to be summed over the quark flavours q, and the subscript i labels up-type quarks 
{i — 1) and down- type quarks (i = 2). The model-dependent coefficients a^i are given by 



giriqi 
AmwB, 



Re {5u [gZ^2 - g'Z^i] )DiCi( \- + 



+ Re {S2i[gZ^2 - g'Z^,]) [ + 



^H, J J 



(3) 
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where 

= Vu^ ^ ^126*^^X1' 

with UiiT^i denoting hypercharge and isospin, and 

Z-^^{Z-^i), S2i ^ Z^4,{-Z-^3) (5) 

5j = sin/3(cos/3), Cj = sin a (cos a), L>j = cosq;(— sina), (6) 

for up (down) type quarks. We denote by mj^j < n^'Hi the two scalar Higgs masses, and a 
denotes the Higgs mixing angle. Finally, we note that the factors rjij arise from the diagonal- 
ization of the squark mass matrices: diagi(m^,m|) = rjM^r)"^, which can be parameterized 
for each flavour / by an angle 9f and phase 7/: 

/ cos 6'/ sin 61/6'^/ \ ^ / r/n 7712 \ 

[- sine fe-'-^f cos0f ) " [ rj2i V22 ) ' ^ ' 

In the models we study below, the squark flavours are diagonalized in the same basis as the 
quarks. 

2.2 Hadronic Matrix Elements 

The scalar part of the cross section can be written as 



4m? 



as = ^ [Zf, + {A- Z)Uf , (8) 

TT 

where is the reduced LSP mass, 

^= E /?,^— + ^/SE— > (9) 

~, nia 27 "Vl nia 

y q=u,a,s y c,b,t y 

the parameters f^^ are defined by 

iT^pfrq = {P\mqqq\p) = rUgBg, (10) 
/tg ~ -'- ~ J2q=u,d,s frq [^^]-> fn has a similar expression. 
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We take the ratios of the quark masses from [19]: 



in TTi 

— = 0.553 ± 0.043, — = 18.9 ±0.8, (11) 
rud rud 



and following [20], we have: 

-B,, — B 



Defining 



we then have 



Bd + 

Bd ^ 2 + {z- l)y 
S„ 2z-{z-l)y' 
The coefficients /t„ are then easily obtained; 



z = i = 1.49. (12) 



y - (13) 



(14) 



_ niuBu _ 2S 



B,i, 



_mdBd _ 2S 



^ 2(^)Ev 
■'^^ m„ m„(l + ^)' ^ ^ 

The final task is to determine the quantity y characterizing the density of ss in the nucleon. 
This may be determined from the 7r-nucleon E term, which is given by 

cr^jv = E = ^(m„ + md)(S„ + 5d)- (18) 

We are motivated to reconsider the value of y in hght of recent re-evaluations of the 7r-nucleon 
sigma term E, which is related to the strange scalar density in the nucleon by 

y = 1 - (7o/E, (19) 

where is the change in the nucleon mass due to the non-zero m, d quark masses, which 
is estimated on the basis of octet baryon mass differences to be do = 36 ± 7 MeV [21]. In 
our previous work [4,5], wc assumed a relatively conservative value S = 45 MeV, which was 
already somewhat larger than naive quark model estimates, and corresponded to ?/ ^ 0.2. 
However, recent determinations of the 7r-nucleon E term have found the following values at 
the Cheng-Dashen point t = +2m^ [13]: 

EcD = (88 ± 15, 71 ± 9, 79 ± 7, 85 ± 5) MeV. (20) 
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These should be corrected by an amount — A^^ — Ao- ~ — 15 MeV to obtain E. Assuming for 
definiteness the value Scd = 79 ± 7 MeV, we may estimate 

S = (64 ±8) MeV. (21) 

Such a relatively large value of S has recently received support from an unexpected quarter, 
namely the apparent observation of exotic baryons 0^,H in an antidecuplet of flavour 
SU(3) [14]. The existence of such states has been a long-standing prediction of chiral-soliton 
models, but the details of their spectroscopy depend, in particular, on the value of E: 

TTl 

— E = 3{4M^ - 3M^ - M^) + 4:{Mn - Ma) 
m , ' , ' 

octet decuplet 
-4(Mh— -Me+) 

antidecuplet (22) 

in the chiral-soliton model. Inserting the nominal values Mq+ = 1540 MeV and M= — = 
1862 MeV, we find S = 72 MeV, corresponding to y ~ 0.5. This determination should be 
taken with a couple of grains of salt, since it is unclear whether either the 0+ or the H 
exist. However, since this value is consistent with the more direct estimate ()21|). we adopt 
S = 64 MeV and 45 MeV as alternative nominal values, corresponding to ?/ ^ 0.45 and 0.2, 
respectively, which we use later to discuss the implications of varying S. 

3 Exploration of the CMSSM 

We begin by considering the constrained version of the MSSM (CMSSM) [22-24]. This 
class of models is fully described by four parameters and a sign: a unified gaugino mass, 
mi/2, a unified scalar mass, mo, a unified trilinear mass term, Aq, and the ratio of the 
Higgs v.e.v.'s, tan/5. In addition, the sign of the /i parameter must also be specified. The 
phenomenology of these models has been well studied. The parameters of models with an 
acceptable cosmological relic density generally fall into one of the following regions: the 
coannihilation region ^, where the mass of the neutralino is nearly degenerate with the mass 
of the stau; the rapid-annihilation funnel, where the mass of the neutralino is close to one- 
half the mass of the pseudoscalar Higgs boson A; and the focus-point region, which is found 

^Note that, because the relic density has now been determined with high accuracy by cosmological 
observations [25], and accelerator limits disfavour low mx/2j we no longer distinguish a bulk region at low 
™i/2 s-nd Too from the coannihilation region. 
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at extremely high values of mo, and is at the edge of the parameter space which allows for 
radiative electroweak symmetry breaking. We start with an examination of some specific 
benchmark parameter choices [6, 17] that populate these allowed regions. 

3.1 Benchmark Scenarios 

Fig. n shows the effect of the value of S on the magnitude of the spin-independent elastic 
X-proton scattering cross section in the specific cases of the CMSSM benchmark scenarios 
discussed in [17]. Points A, B, C, D, G, H, I, J, and L are in the coannihilation region, points 
K and M are in the rapid-annihilation funnels, and points E and F are in the focus-point 
region. Point M is not shown as its cross section falls below the scale of the plot ^. It is 
clear that the value of S has quite significant impact in all the scenarios, as indicated by the 
behaviours of the different lines. There is a general trend for the cross section to increase 
approximately quadratically with S. This would be exact if the < p\uu\p > and < p\dd\p > 
contributions were negligible compared with the < > contribution. However, Fig. ^ 

shows that the increasing trend is not exactly universal, reflecting the different relative 
weights of the various < p\qq\p > contributions in the different benchmark scenarios. These 
depend on tan/5 and the sign of the Higgs-mixing parameter /i, as can be seen from the 
formulae in the previous Section. 

We have plotted in Fig. Rvalues of the cross section corresponding to S > 36 MeV, i.e., 
consistent with assuming y > 0. The dashed curves in Fig. ^ around benchmark points C, 
G, and J display the effect of the uncertainty in cxo as well as the mass ratios which enter 
into the determination of the /j-, and ultimately the elastic cross section. We see that this 
uncertainty is not negligible, although that associated with S is clearly more important. 

We see that, in all scenarios and for all plausible values of S, the estimated cross section 
lies considerably below the current upper limits of CDMS II [10], which can at best exclude 
models with cross sections larger than 3 x lO^'' pb when = 60 GeV. If future experiments 
achieve a sensitivity of 10^^ pb, one can plainly see that several of the benchmark scenarios 
will be probed, particularly if S is large. 

It is clear from the above discussion that better understanding of the non-perturbative 
hadronic matrix elements S and ctq will be needed before the spin-independent elastic- 
scattering cross section can be predicted accurately in any specific supersymmetric model. 
This means, in particular, that unless these hadronic matrix elements can be determined more 
accurately, it will be difficult to convert any LHC or LC measurements of MSSM parameters 

^ These benchmark points were formulated assuming mt = 175 GeV. The small shifts required if one uses 
the new central value mt — 178 GeV do not impact significantly the cross sections calculated in Fig. 
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Figure 1: The dependences on the ir-nucleon E term of the elastic cross sections of the 
benchmark points [17]. The dashed lines indicate the sensitivities to ctq i'^ the cases of 
benchmark scenarios C, G and J. The predicted cross sections are smaller than the CDMS II 
upper limits [10] for the the models considered, for all E values shown. 

into accurate predictions for elastic-scattering rates, even if they do suffice to calculate accu- 
rately the relic LSP density. The experimental determination of S is notoriously uncertain: 
perhaps the time is ripe for another lattice QCD approach? 

The benchmark scenarios discussed above were formulated within the CMSSM, and our 
next step is to explore the CMSSM more generally, to see whether larger cross sections are 
possible in regions of its parameter space. 

3.2 General Analysis of CMSSM Models compatible with WMAP 

As is well known, for any given value of tan/3, Aq and rrit, the CMSSM parameter space 
consists of narrow strips in the (mi/2,mo) plane, where the relic density falls within the 
range allowed for cold dark matter by WMAP and other experiments. In the following, we 
no longer consider results in the focus-point region: this now appears at very large mo if 
one adopts the new central value rrit = 178 GeV ^, as we do henceforth. At low values of 
mi/2, the length of the strip is in turn restricted by experimental constraints such as rrih, 
m^± and b s'-f, whereas at high values of mi/2 the strips are truncated by the relic density. 
We display in Fig. |21 the {rni/2,rno) planes for tan/3 = 10 and (a) /i < 0, (b) /i > 0, (c) 
tan/3 = 40, /i < and (d) tan/3 = 57, fi > 0. The latter choices of tan/? are close to the 
•^We use ■mh{mb)jjg — 4.25 GeV throughout. 
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maximal values we now find for the corresponding signs of /i. These have increased with the 
new best-fit value rrit = 178 GeV and recent improvements in our spectrum evaluation code ^. 
The rapid-annihilation funnels visible in panels (c, d) are located at values of mi/2 that are 
similar to what we would have found previously for tan f3 = 35, /i < and tan (3 = 50, /i > 0. 

The various experimental and cosmological constraints on the CMSSM are displayed in 
various {mi/2,mQ) planes in Fig. El but we do not use them all as absolute hmits. The dark, 
tan-shaded regions are, however, completely excluded because there the LSP is charged, 
being the lighter f. The thin blue strips are those favoured by the WMAP constraint on 
the rehc density of cold dark matter: 0.094 < VtcDuh^ < 0.125 if VL^ ~ ^cdm, and we also 
display the restrictions that the accelerator constraints due to uih (red dash-dotted lines), 
m^± (black dashed lines) and b ^ s'j (medium, green shading) impose on the ranges of mi/2 
and hence mg allowed along the WMAP strips. The constraints that would be imposed by 
(7^ — 2 at the 2-cr level if the Standard Model contribution is evaluated using e'^e~ annihilation 
data alone, neglecting r decay data, are shown by light, pink shading in panels (b, d) ^. 

Each of the panels also displays contours of the spin-independent elastic-scattering cross 
section calculated for S = 45 (lighter, blue dashed contours) and 64 MeV (black dotted con- 
tours), labelled by their exponents in units of picobarns. We see that, for /i > in panels (b, 
d) of Fig. the cross-section contours progress monotonically downward as mi/2 increases, 
with the E = 45 MeV contours always at smaller mi/2 than the corresponding contours for 
S = 64 MeV. However, the progression is not monotonic for /i < 0, as seen in panels (a, 
c). This is because of the possibility of a cancellation between different contributions to the 
scattering amplitude [4]. 

For the purpose of this paper, we choose to treat the WMAP constraint as an upper 
limit on fi^^/i^ = fx^CDAih"^ : /x ^ 1; thus allowing for another component of cold dark 
matter with a fractional density 1 — > 0. In this case, the small regions of the (mi/2, '^o) 
planes between the WMAP strips and the charged LSP corners are also allowed. We can 
see in Fig. |21 that the spin-independent elastic scattering cross section is very similar in the 
underdense regions with f-^ < 1, which lie below the WMAP strips and above the charged 
dark matter region, and also those inside the rapid-annihilation funnels for large tan (3. 

Implementing the accelerator constraints, using the relic density allowed by WMAP as an 

upper limit: fl^h"^ = fx^cDAdh"^, and rescaling the cross section by a factor if fx < 1, so 

as to account for the fact that neutralinos could constitute only a fraction of the galactic 

^The most recent improvements include implementation of the full set of two-loop rcnormalization-group 
equations. 

'"^See [16] for a discussion on the g^~2 deviation range used here. We recall that no models with yu < 
would be allowed at this significance level. 
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Figure 2: The (mi/2,mo) planes in the CMSSM for (a) tan/5 =10, /i < 0, (b) tan/3 = 
10, /i > 0, (c) tan/3 = 40, fi < and (d) tan/3 = 57, /i > 0, a// assuming Aq = 0. VFe 
display the WMAP relic-density constraint, the experimental constraints due to mh, ni-^±, 
b ^ s'j and g^^ — 2, and contours of the spin-independent elastic-scattering cross section 
calculated for S = 45 and 64 MeV (lighter, blue and black dotted contours, respectively), 
labelled by their exponents in units of picobarns. 
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halo and that there would be another important local component of cold dark matter, we 
find the ranges for the effective spin-independent elastic-scattering cross section shown in 
Fig. El These ranges were obtained by statistical sampling of the allowed regions of the 
CMSSM parameter spaces for the indicated parameter values. The sampling was performed 
over values of mi/2 = 0.1 to 2 TeV, mo = to 2 TeV, tan/5 = 2 to 43 (58) for fi < (>)0, 
and Aq = —3 to +3mi/2- Because of the rescaling and the fact that regions with < 1 
have similar intrinsic cross sections to regions with = 1, the points with in the range 
of ^cDM favoured by WMAP generally appear at the top of the allowed ranges. In general, 
the calculated cross sections lie below the present CDMS II upper limits, except for certain 
models with the smallest values of mi/2 that are allowed when tan/3 ~ 10 and > 0, if one 
uses S = 64 MeV. 

3.3 Preferred Range of Sparticle Masses 

Progressing beyond the above implementation of laboratory experimental constraints, the 

sparticle mass range preferred within the CMSSM has recently been reassessed [16], in light 

of recent precision measurements and higher-order calculations in the Standard Model and 

the MSSM. As has already been recalled, the anomalous magnetic moment of the muon, 

— 2, disagrees with the Standard Model by between 2.5 and 3 standard deviations [26], 

if low-energy e+e~ data are used to estimate the strong-interaction contribution to — 2: 

see the pink shaded regions in Fig. El The central experimental value favours fi > and 

mi/2 ~ 300 GeV for tan/3 = 10 ^, and the preferred value of mi/2 increases with tan/3. 

The present central values of Myy and sin^ O^f / also disagree marginally with the latest 

theoretical calculations within the Standard Model. Given the errors, these discrepancies 

are not significant in themselves, but it so happens that they are each, separately, best fit 

also by mi/2 ~ 300 GeV for tan/3 = 10. The quality of fit in the {mi/2,A) planes for 

tan /3 = 10, 50 has been explored, and the 68% and 90% confidence level regions have been 

delineated: they stretch up to mi/2 ^ 1000 GeV [16]. 

In Fig. m we display scatter plots of the spin- independent elastic-scattering cross section 

calculated for S = 45 and 64 MeV, as usual, for the portions of the WMAP strips allowed 

for (a, b) tan /3 = 10, /i > and (c, d) tan /3 = 50, /i > at both the 68% and 90% confidence 

levels. The two choices CL = 68 % and 90 % have different colours (darker, blue x and 

lighter, green + signs, respectively). We do not see large qualitative differences between the 

cross-section predictions in the 68% and 90% confidence-level cases. Also, comparing the 

^We note in passing that the minimum of the function almost coincides with benchmark point B 
of [17], to which Point la of [27] is also similar. 
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Figure 3: Scatter plots of the spin-independent elastic- scattering cross section predicted in 
the CMSSM for (a, b) < and (c, d) fi > 0, with (a, c) T, = 45 MeV and (b, d) 
S = 64 MeV. 
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top and bottom panels, we do not see large qualitative differences between the cross-section 
predictions in the cases tan (3 = 10, 50, though the latter are slightly larger. However, 
comparing the left and right panels, we once again see the direct effect on the cross section 
due to our choice of S. Since /z > in this analysis, there is no possibility of a cancellation in 
the cross section. Moreover, comparing with the corresponding panels of Fig. El we note that 
the preferred range of mi/2 and hence happens to be that where the spin-independent 
elastic scattering cross section is relatively large. 

We see that an improvement in the present CDMS II limit by an order of magnitude 
would just begin to touch the estimated cross-section range, for low and large S. On the 
other hand, an improvement by around 4 orders of magnitude would be required to cover 
completely all the regions allowed at the 90% confidence level for all the considered range 
45 MeV< S <64 MeV. 

4 Detection in Models with Non-Universal Scalar Masses 

Larger cross sections may be found in models in which the CMSSM assumptions of universal 
soft supersymmetry-breaking scalar masses ttlq are relaxed, as we now discuss. 

4.1 General Low-Energy Effective Supersymmetric Theory 

We first consider relaxing the universality assumption for the Higgs bosons and for the soft 
supersymmetry-breaking squark masses relative to those of the sleptons, requiring only that 
all the squark and slepton squared masses remain positive under renormalization up to the 
GUT scale. This we term the most general low-energy effective supersymmetric theory 
(LEEST) [28]. It is clear that relaxing the CMSSM relationship between the squark and 
slepton masses rrig, might have a direct impact on the elastic-scattering cross section, 
although the freedom to adjust niq/mi is severely restricted by the LEEST requirement that 
the squared masses remain positive up to the GUT scale. The primary impact of relaxing 
universality for the Higgs boson masses is to permit variations from the CMSSM values of 
the pseudoscalar Higgs mass and the magnitude of Higgs mixing |yu|, which are fixed by 
the electroweak vacuum conditions. We discuss below the extent to which these different 
effects can be disentangled. 

We display in Fig. El scatter plots of the spin- independent elastic-scattering cross section 
for both signs of the Higgs-mixing parameter /i: negative in panels (a, b) and positive in 
panels (c, d). Predictions for two values of S, the conservative value of 45 MeV and the 
more modern value of 64 MeV, are shown in panels (a, c) and (b, d), respectively. We see 
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Figure 4: Scatter plots of the spin-independent elastic-scattering cross section predicted in the 
CMSSM for (a, h) tan/? = 10,/i > and (c, d) tan/? = 50,/i > 0, with (a, c) J: = A5 MeV 
and (h, S = 64 MeV. The predictions for models allowed at the 68% (90%) confidence 
levels are shown by blue x signs (green + signs). 
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that predictions of S for /i < never rise to the sensitivity of the CDMS II experiment [10], 
whichever value of S is used. However, a few points with ^ 700 GeV do exceed the 
current CDMS II limit for /i > 0, as seen in panels (c, d), particularly when the larger value 
of S is used. We discuss the nature of these excluded points further below. 

4.2 Models with Non-Universal Higgs Masses 

Since the parameter space of the LEEST has quite a large dimensionality, it is difficult 
to visualize clearly what classes of models might be excluded by CDMS II. This becomes 
clearer if one considers a class of models with a lower-dimensional parameter space, namely 
those with universal soft supersymmetry-breaking masses for squarks and sleptons but non- 
universal Higgs masses (NUHM) [5,29,30], which allow values of and ttla differing from 
those in the CMSSM. We display in Fig. El scatter plots of the spin- independent elastic- 
scattering cross section for (a, b) /i < and (c, d) /i > 0. The same two choices of S, namely 
45 MeV and 64 MeV are made in panels (a, c) and (b, d), respectively. 

We see again that no NUHM points can be excluded for fi < but that, as in the LEEST 
case, some /i > NUHM points may be excluded by CDMS II. This is true, in particular, for 
the larger choice of S. The similarities between the general trends in the corresponding panels 
of Figs, m and El suggest that the dominant effects may be due to relaxing the universality 
assumption for the Higgs masses which, we recall, allows the values of and to differ 
from those in the CMSSM. In fact, the LEEST does not have much leeway for varying the 
ratio mq/mi at low energies, since we restrict the soft supersymmetry-breaking scalar masses 
so that the effective scalar squared masses remain non-tachyonic all the way up the GUT 
scale. 

The fact that the most significant variations from the CMSSM are likely to be those in 
and rriA is supported by a previous general study of the NUHM [30], in which various 
(mi/2, '^o)) {^A, fJ') and (m^, M2) planes were exhibited. The behaviours of the cross section 
in the NUHM (mi/2,''^o) planes were similar to those found in the CMSSM, varying mainly 
with mi/2 and less with niQ [5]. The dependence on M2 in the (mA, M2) planes basically 
refiected the same mi/2 dependence. The most striking dependence of the cross section was 
on so we focus here on the (m^,//) planes for mi/2 = 500 GeV, niQ = 1000 GeV and 
different choices of tan (3, which are displayed in Fig. 13 Regions outside and below the 
black double-dash-dotted lines have negative Higgs masses-squared below the GUT scale, 
and are hence unstable, so only the regions between and above these lines are allowed. This 
constraint becomes less important as tan j3 is increased. 
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Figure 7: Contours of the spin-independent elastic cross section in the {mA, fi) planes for 
(mi/2,mo) = (500, 1000) GeVandtanp = (a) 10, (h) 20, (c) 35 and (d) 50, for^ = 45 MeV 
(dashed blue lines) and S = 64 MeV (black dotted lines), labelled by their exponents in units 
of picobarns. The regions excluded by CDMS II [10] lie below the solid black lines. 
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As usual, the dark, green shaded regions are excluded by 6 — > 57, the light, pink shaded 
regions are those preferred by (7^ — 2, the solid dark blue strips are those where falls 
within the range preferred by WMAP, the red dash-dotted line is the Higgs mass constraint 
and the black dashed line is that imposed by the chargino mass. The outward bulges in the 
WMAP strips are caused by rapid-annihilation funnels. The Higgs constraint forbids regions 
with low mA, which are also excluded by the GUT Higgs stability constraint for tan f3 = 10, 
as seen in Fig. [7|^a), but not necessarily for larger values of tan /?, as seen in the other panels 
of Fig. [7| The chargino constraint removes regions with small |yu|. 

Contours of the spin- independent elastic cross section are also plotted in the (m^,/i) 
planes for various values of tan j3 in Fig. El labelled by the exponents in units of picobarns 
(blue dashed curves for S = 45 MeV, black dotted curves for S = 64 MeV). We see that the 
largest values of the spin-independent elastic scattering cross section occur when /z and 
are relatively small. Also displayed in Fig. [3 are the regions excluded by the CDMS II upper 
limit (solid black line), including also the factor < 1 where appropriate for models with 

< VLcdm- In panel (a) for tan/3 = 10, the regions excluded by CDMS II were already 
excluded by the GUT Higgs stability and 6 — > 37 constraints. However, in the other panels 
we see that there are regions at low /i and tua that were allowed by the other constraints 
but are now excluded by CDMS II. These regions become progressively more extensive as 
tan (3 increases. 

These regions are reflected in Fig. |Hfa), which displays in the {fi,mA) plane the NUHM 
points from Fig. IH^c, d) that are excluded by the CDMS II result if one assumes S = 64 MeV 
(dark, red x signs) or S = 45 MeV (lighter, blue squares). As expected, they cluster at 
small values of fi and tua ^. Their values of /i and are generally smaller than those of 
the benchmark points [17], which are all compatible with CDMS II, as we saw in Fig. d For 
comparison, only benchmark point B (m^ ~ 370 GeV) has a pseudoscalar mass less than 
400 GeV and all but points B, I (m^ ~ 450 GeV), and L {rriA ~ 490 GeV) have pseudoscalar 
masses in excess of 500 GeV. Similarly, with the exception of the focus points (E and F), 
typical values /i are relatively large. Point B has ^ ~ 330 GeV, point I has ~ 440 GeV, 
and G has /i ~ 470 GeV, whereas all other points have fi in excess of 500 GeV. Fig. |Hl^b) is 
the corresponding plot for the excluded LEEST points from Fig. El^c, d). This exhibits very 
similar features, confirming the importance of these variables also in the LEEST scenario. 
In contrast, the ratios m^/m^ for the excluded LEEST points do not exhibit any clustering 
at low values. 

^Analogous high-cross-section points for /i < are excluded by the 6 — > 57 constraint, as seen in Fig. [3 
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Figure 8: Scatter plots in the {fi,mA) plane of points from (a) Fig.\^c, d) and (h) Fig.\^c, 
d) that are excluded by the CDMS II constraint for S = 64 MeV (dark, red x signs) or 
S = 45 MeV (lighter, blue squares), exhibiting similar clustering at low values of fi and mA- 
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5 Conclusions and Prospects 



In this paper, we have made a new comparison between theoretical predictions of the spin- 
independent cross section for the elastic scattering of supersymmetric dark matter and the 
improved experimental upper limit recently provided by CDMS II [10]. In making this 
comparison, we have contrasted the theoretical predictions made with different estimates of 
the vr-nucleon S term. Larger values may be supported by recent reports of exotic baryons, 
but these do not increase greatly the ranges of theoretical models excluded by CDMS II. We 
have also incorporated in our analysis the new central value of rrit, which enters indirectly 
into constraints on the supersymmetric parameter space and into relic-density calculations. 

Some supersymmetric models with non-universal Higgs masses (NUHM) are now ex- 
cluded by the CDMS II upper limit, as are some models which also incorporate non- universal 
squark and slepton masses (LEEST). These are mainly models with the smaller values of 
and/or that become allowed when the universality conditions are relaxed for the Higgs 
masses. 

On the other hand, only very small parts of the CMSSM parameter space are yet excluded. 
Specifically, the cross sections we find in the supersymmetric benchmark scenarios of [17] all 
lie considerably below the CDMS II sensitivity, as do all points allowed at the 68% or even 
90% confidence level by a recent likelihood [16] analysis of the CMSSM parameter space 
incorporating information on niw, sin^ 9eff and — 

An improvement over the present CDMS II sensitivity by about an order of magnitude 
would begin to challenge the preferred region of CMSSM parameter space, but an improve- 
ment by about four orders of magnitude would be required to cover it completely. We 
conclude that direct searches for supersymmetric dark matter are just beginning to reach 
interesting sensitivities, but that considerable improvement will be needed to exclude (or 
hopefully discover) supersymmetric dark matter. 
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